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PREFACE 


This  is  one  of  a  series  of  reports  that  summarize  data  relating  to 
blasting  procedures  and  blast  effects  in  cold  regions.  They  are  organized 
to  deal  successively  with  (1)  explosions  in  air,  (2)  explosions  in  water, 

(3)  explosions  in  solid  ground  materials.  For  the  most  part,  the  blasting 
procedures  used  in  cold  regions  are  not  much  different  from  those  that  are 
in  general  use  elsewhere.  Similarly,  the  principles  involved  in  assessing 
blast  effects  in  cold  regions  are  the  same  as  those  that  apply  to  blast  ef¬ 
fects  generally.  The  reports  therefore  summarize  principles  and  data  for 
general  explosions  technology,  and  then  present  the  procedures  and  data  for 
cold  environments  within  this  framework. 

Tlie  purpose  of  the  series  Is  to  provide  a  convenient  reference  source 
for  engineers  faced  with  problerus  of  explosions  or  blasting  in  cold  re¬ 
gions.  Secause  not  all  engineers  arc  well  acquainted  with  explosions  tech¬ 
nology.  relevant  physical  principles  are  explained  or  summarized,  but  there 
is  no  attempt  to  explore  the  underlying  theory  in  depth,  not  is  there  any 
treatment  of  the  practical  aspects  of  explosives  use  and  safety  practices. 
These  topics  are  covered  well  in  Army  Technical  Manuals  and  Army  Materiel 
Command  publications,  as  well  as  in  commercial  blasters'  handbooks  and  In 
textbooks. 

This  report  was  prepared  by  Dr.  Malcolm  Mellor,  Research  Physical  Sci¬ 
entist,  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions  Research 
and  Engineering  Laboratory.  The  work  was  done  under  DA  Project  4A762730- 
AT42,  Design,  Construction,  and  Operations  Technology  for  Cold  Regions, 

Task  Area  Combat  Support,  Work  Unit  029,  Explosives  and  Projectile  Inpact 
Under  Winter  Conditions. 

The  author  Is  grateful  to  Dr«  Llndamae  Peck,  and  Donald  Mbert  for 
their  careful  revxews  of  the  manuscript  and  suggestions  for  improvement. 

The  numerous  figures  were  prepared  by  Edward  Perkins  and  his  staff. 
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BLASTING  AND  BLAST  EFFKCTS  IN  COLD  REGIONS 
PART  I:  AIR  BLAST 

Malcolm  Mellor 


Introduction 

When  a  shock  wave  propagates  directly  through  free  air,  there  is  not 
much  dift=»rence  between  air  blast  in  cold  regions  and  air  blast  anywhere 
else.  Shock  properties  are  influenced  by  differences  of  temperature,  air 
density,  ind  atmospheric  pressure,  but  even  in  temperate  regions  there  are 
large  variations  of  temperature,  pressure  and  density  with  altitude,  sea¬ 
son,  and  geographical  location.  In  short,  the  procedures  that  exist  for 
assessing  the  effects  of  varying  atmospheric  conditions  in  ten^erate  re¬ 
gions  are  fully  adequate  for  dealing  with  the  conditicns  encountered  at 
high  latitudes  and  high  surface  elevations, 

A  special  feature  of  cold  regions  is  widespread  occurrence  of  snow 
cover,  ice-covered  water  bodies,  hard  glacial  ice,  and  frozen  ground  sur¬ 
faces.  Of  these,  snow,  and  possibly  floating  ice,  can  be  expected  to  re¬ 
flect  strong  air  shocks  differently  than  hard  surfaces,  with  consequent  ef 
fects  on  the  spatial  distribution  of  shock  properties. 

The  characteristic  ground  surfaces  of  cold  regions  can  have  signifi¬ 
cantly  different  thermal  and  optical  properties  from  those  of  common  sur¬ 
faces  in  temperate  regions.  This  has  potential  significance  with  respect 
to  air  blast  effects  from  nuclear  explosions,  but  it  seems  that  snow,  ice 
and  frozen  ground  tend  to  limit  the  effects  of  thermal  radiation,  making  a 
nuclear  blast  behave  more  like  a  conventional  explosion,  vn  ,h  formation  of 
a  precursor  unlikely.  In  other  words,  these  surfaces  are  "near-ideal," 
like  water,  light-colored  concrete,  and  packed  soil  with  sparse  vegetation 

Falling  snow,  blowing  snow,  and  ice  fog  can  attenuate  optical  and 
thermal  radiation  from  a  nuclear  explosion,  as  can  rain  and  fog  in  warmer 
conditions.  The  optical  properties  of  all  these  suspensions  are  known  in 
sufficient  detail  for  calcul  ion  of  the  transmission  and  attenuation  of 
radiation,  but  this  topic  is  outside  the  scope  of  the  present  report. 


All  of  uiiLr?  suggests  that  the  main  requirement  for  understanding  air 
blast  In  cold  regions  Is  a  knowledge  of  standard  air  blast  technology,  to¬ 
gether  with  some  specific  Information  oi  flection  and  attenuation  over 
snow  surfaces. 

Ideal  Blast  Waves  in  Free  Air 

ivn  e?;ploslon  In  air  forms  gas  at  very  high  pressure  and  compresses  the 
surrounding  air,  leading  to  expansion  at  supersonic  speeds  and  to  conse¬ 
quent  formation  of  a  shock  front.  If  the  air  surrounding  the  explosion  is 
homogeneous  and  of  unlimited  extent,  the  gas  bubble  expands  adlabatlcally 
in  a  symmetric  spherical  pattern.  Pressure  at  the  shock  front  rises  from 
ambient  to  a  high  value,  called  the  peak  overpressure ,  almost  Instantane¬ 
ously.  Pressure  behind  the  shock  front  decreases  as  the  gas  bubble  expands 
and  changes  temperature.  It  can  drop  below  ambient  atmospheric  pressure 
because  of  Inertia,  and  so  pulsate  to  give  secondary  and  tertiary  pulses. 

As  the  spherical  shock  front  expands  and  the  pressure  decreases,  the  shock 
velocity  decreases  until  the  front  eventually  becomes  an  acoustic  wave. 

For  an  observer,  or  a  gauge,  at  a  fixed  distance  from  the  explosion, 
pressure  varies  with  time  as  shown  In  Figure  1.  Between  the  instant  of 
detonation  and  the  arrival  of  the  shock  front,  pressure  remains  at  the  am¬ 
bient  level.  It  then  increases  to  a  peak  value  almost  Instantaneously,  and 
as  the  shock  front  passes  the  pressure  decreases  with  time.  Looking  at 


Figure  1.  Variation  of  overpressure  with 
time  for  an  air  shock  passing  a  fixed  ob¬ 
servation  point. 
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things  another  way,  as  In  Figure  2,  the  pressure  at  any  given  Instant  Is  a 
maxlraura  at  the  position  of  the  shock  front,  and  it  decreases  In  the  direc¬ 
tion  of  the  source  oi  the  explosion.  If  these  pressure  curves  are  plotted 
for  successive  instants,  the  locus  of  the  peak  pressure  la  the  attenuation 
curve,  i.e.  the  decay  of  peak  overpressure  with  distance  from  the  source. 
The  jump  in  pressure  across  the  shock  front  is  accompanied  by  jumps  in  den 
slty,  particle  velocity  and  temperature.  The  variations  of  density  and 
particle  velocity  with  distance  and  time  are  quali-tatlvely  similar  to  the 
variation  of  overpressure  as  shown  in  Figure  2. 

If  the  explosion  is  one  that  produces  a  relatively  slow  pressure  in¬ 
crease,  say  from  a  gas  discharge  or  from  the  burning  of  low  explosive,  the 
wave  front  near  the  source  may  have  a  'elatlvely  long  rise  time.  However, 


0i«fQnc6  from  Caplosion 


Figure  2.  Variation  of  overpressure  with  distance  at 
a  given  instant.  Diagrams  are  drawn  for  several  dif¬ 
ferent  times,  so  that  the  locus  of  peak  values  gives 
the  attenuation  of  peak  overpressure  with  distance. 


Direction  ot  PrgpogotiOn 


Distance 


Figure  3.  Change  of  wave  form  as  a  pressure  wave  with  long 
initial  rise-time  steepens  with  time  and  distance  so  as  to 
transform  into  a  true  shock.  Attenuation  is  ignored  in  this 
diagram. 


the  front  can  steepen  as  It  moves  out,  becoming  a  true  shock  at  some  dis¬ 
tance  from  the  source  (Fig.  3).  This  effect  Is  sometimes  referred  to  as 
"shocking  up." 

Overpressure  is  the  difference  between  ambient  pressure  and  pressure 
behind  the  shock  front.  There  Is  also  a  sort  of  wind  pressure,  called  the 
dynamic  pressure ,  associated  with  passage  of  the  blast  wave.  It  is  deter¬ 
mined  by  air  density  p  behind  the  shock  front  and  by  the  particle  velocity 
u; 

q  =  ^  0  u^  (  1) 


where  q  Is  the  dynamic  pressure.  When  a  solid  obstacle  is  struck  by  air 
blast,  the  Incident  and  reflected  overpressures  of  the  shock  front  at  first 
dominate  the  loading  process,  producing  diffraction  effects  and  creating 
vortices  something  like  those  formed  by  strong  winds.  After  passage  of  the 
shock  front  the  dynamic  pressure  loads  the  obstacle  like  a  wind.  This 
"wind"  can  reverse  when  the  positive  phase  passes  through  and  Is  followed 
by  the  negative  phase  of  the  blast  wave. 

The  Shock  Equations  for  Air  Blast 

Shock  equations  can  bo  derived  for  the  process  of  detonation  In  aii  ex¬ 
plosive,  for  the  compression  and  displacement  of  air  around  an  explosion 
source,  or  for  the  mechanical  response  of  any  liquid  or  solid  medium  sur¬ 
rounding  an  explosion.  The  same  theory  is  applied  in  all  th’se  cases  to 
develop  the  shock  relations,  or  Rankine-Hugoniot  conditions.  By  consider¬ 
ing  the  conservation  of  mass,  momentum  and  energy  across  the  "jump"  repre¬ 
sented  by  the  shock  front,  and  taking  a  coordinate  system  that  travels  at 
the  speed  of  the  wave  front  (Lagrange),  the  following  relations  are  ob¬ 
tained  : 


Conservation  of  mass 
Conservation  of  momentum 

Conservation  of  energy 


U2P2  =  UiPi  (2) 

2  2 

P2  +  P2U2  =  Pi  +  PiUi  (3) 
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where  p,  p,  and  u  are  pressure,  density  and  particle  velocity  respectively, 
and  e  is  the  internal  energy.  The  subscripts  1  and  2  refer  to  conditions 
on  each  side  of  the  shock  front.  For  a  stationary  (Euler)  coordinate  sys- 
tetii,  In  which  the  shock  front  travels  at  velocity  U,  the  equations  become: 
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p  +  pCu-u)*"  ”  Po  +  Po(uo“^)^ 

P  (|-  u^  +  e)(u-u)  +  pu  =  Po(^  “0^  +  eo)(uo-U)  +  pouQ-  (7) 

The  subscript  "0"  refers  to  ambient  conditions  before  the  arrival  of  the 
shock,  and  the  symbols  without  subscripts  apply  Immediately  behind  the 
shock.  In  most  cases  the  ambient  air  Is  stationary,  so  that  uq  la  zero. 

To  these  mechanical  conditions  are  added  the  Ideal  gas  laws: 


pV 


p/p 


RT 


(8) 


where  R  is  the  gas  constant,  T  Is  absolute  temperature  and  V  is  specific 
volume.  For  a  polytropic  gas  the  Internal  energy  e  is 


c  T  =  - -Y  pV 

V  Y~f 


_i_  £ 

y-l  p 


(9) 


where  Cy  Is  the  specific  heat  at  constant  volume  and  y  is  the  adiabatic 
exponent,  equal  to  the  ratio  of  specific  heat  at  constant  pressure  to  spe¬ 
cific  heat  at  constant  volume,  Cp/cy.  For  adiabatic  compression  or  ex¬ 
pansion 


pV 


=  E_ 


A  = 


RT 


(10) 


P  P' 

where  A  Is  a  constant.  The  speed  of  sound  In  the  gas  els 


2  _  dp 

"  dp  ~ 


yRT  =  ypV 


YP/P 


(11) 


Combining  these  relations,  useful  expressions  for  various  properties 
of  the  shock  front  can  be  obtained.  These  expressions  are  given  later. 

Scaling  Procedures  for  Comparison  of  Explosions 

One  of  the  numerous  variables  In  explosions  is  the  size  of  the  charge. 
To  simplify  the  comparison  of  explosions  of  different  size,  it  Is  desirable 
to  find  some  way  to  scale  the  variables  with  respect  to  charge  size. 

For  a  given  type  of  explosive,  pressure  and  detonation  velocity  at  the 
source  are  more  or  less  Invariant  with  charge  size,  so  the  significant  var¬ 
iables  for  the  charge  are  energy  content  and  physical  size.  If  the 
density,  specific  energy  and  shape  of  the  charge  are  constant,  then  charge 
volume  and  energy  are  directly  proportional.  The  simplest  scaling  would  be 
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based  on  geometric  similitude,  and  this  turns  out  to  be  applicable  over  a 
wide  range  when  the  given  medium  surrounding  the  explosion  is  homogeneous. 
Thus  linear  dimensions  of  the  problem  can  be  scaled  with  respect  to  a  line¬ 
ar  dimension  of  the  charge,  to  the  cube  root  of  energy  content,  or  to  the 
cube  root  of  charge  weight.  The  cube  root  version  of  this  scaling  Is  some¬ 
times  referred  to  as  Hopklnson  scaling.  Detailed  and  perceptive  discus¬ 
sions  of  scaling  have  been  given  by  Jones  (1962)  and  Baker  (1973). 

Cube-root,  or  Hopklnson,  scaling  can  be  Illustrated  by  considering  the 
blast  propagation  from  spherical  charges  of  different  size  In  an  infinite 
homogeneous  medium  with  unvarying  properties  (Fig.  4).  If  we  arc  Interest¬ 
ed  In  the  radius  r  at  which  peak  overpressure  has  a  certain  value  p  In  both 
explosions,  then 
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since  the  energy  yield  of  the  charge  E  Is 

4 

E  “  -7  ttR  "  e’ 
3  c 


and  the  charge  weight  is 
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3"\  Pg 


(12) 
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where  e'  is  energy  per  unit  volume  and  pg  Is  unit  weight,  both  constant. 
Rj,  is  the  charge  radius.  This  means  that  measurements  of  peak  overpres- 
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Figure  4.  Scaling  of  shock  wave  properties 
with  respect  to  the  size  and  energy  of  the 
charge 
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sure  p  versus  radius  r  for  exploslotis  of  different  size  will  form  a  single 
curve  when  p  Is  plotted  against  r/R^,,  against  r/E^/3^  oj.  against 
r/wl/3. 

Because  the  shock  velocity  U  Is  related  directly  to  peak  overpressure, 
the  arrival  time  t^  scales  in  the  same  way  as  radius: 
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(15) 


The  positive  phase  duration  tp  (see  Fig.  1)  scales  In  exactly  the  same 
way.  Thus  characteristic  times  for  a  shock  front  with  a  given  pres  •u»‘e  p 
are  scaled  as  t/R^  or  t/El/^  or  t/wl/^. 

Another  quantity  that  scales  In  the  same  way  as  distance  and  time  is 
the  Impulse  1,  which  Is  the  Integral  of  pressure  with  respect  to  time  for 
the  positive  or  negative  phase  of  the  shocli.  For  the  positive  phase  dura¬ 
tion  tp, 

t  +t 
a  p 

I  =  /  (p  -  Po)dt  (16) 

t 

a 

where  p  Is  a  function  of  time,  pg  is  ambient  pressure,  and  t^  is  arrival 
time  for  the  shock  front.  For  scaling. 
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Of  the  various  scaling  factors,  R^.  is  the  most  attractive  in  the  ab¬ 
stract,  since  it  provides  the  basis  for  forming  dimensionless  variables. 
However,  the  cube  root  of  charge  weight  is  far  more  attractive  for  practi¬ 
cal  work,  since  weight  is  often  the  only  thing  the  user  knows  about  the 
charge.  Figure  5  gives  the  conversion  between  "R^,  scaling"  and  "Wl/^ 
scaling. " 

Strictly  speaking,  explosive  type,  specific  enervy  j  d  charge  density 
have  to  be  constants  for  scaling  to  apply,  but  in  prac^-rte  this  is  not  of 
major  importance.  Among  typical  explosives,  speciHc  energy  and  density 
vary  within  a  fairly  narrow  range,  say  by  a  factor  of  N,  where  1.0  <  N  < 
1.5,  For  quantities  scaled  with  respect  to  F^/^,  the  effect  on  the  scal¬ 
ing  factor  is  N^/^,  or  less  than  15%. 
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Figure  5.  Conversion  of  scaled  dimensions. 

The  vertical  axis  Is  length  normalized  with 
respect  to  charge  radius.  The  horizontal 
axis  is  length  scaled  with  respect  to  the 
cube  root  of  charge  weight,  for  both  tradi¬ 
tional  and  SI  units.  The  four  lines  give 
the  conversion  for  four  different  values  of 
charge  density. 

Standard  curves  normally  use  variables  that  are  scaled  with  respect  to 
the  specific  energy  of  TNT.  If  the  data  are  to  be  applied  to  explosive.. 
other  than  TNT,  or  if  data  for  other  explosives  are  to  be  plotted  on  stand¬ 
ard  curves,  a  weight  adjustment  can  be  made.  The  equivalency  factors  dif¬ 
fer  according  to  the  property  considered  and  the  medium  involved.  For  air- 
blast  overpressure,  equivalent  weight  factors  are  given  for  a  range  of  ex¬ 
plosives  in  Table  1.  For  comparison  with  nuclear  explosives,  W  for  the  nu¬ 
clear  explosion  is  often  taken  as  0.5  times  the  radiochemical  yield  ex¬ 
pressed  in  kilotons,  i.e.  1  kiloton  would  be  equivalent  to  10^  lb  of  TNT 
for  generation  of  airblast  overpressure. 

If  the  ambient  pressure  of  the  homogeneous  atmosphere  is  different  for 
different  explosions,  a  different  set  of  scaling  relations  is  used.  In 
this  schemii,  knovm  as  Sachs  scaling,  dimensional  analysis  is  used  to  derive 
normalizing  parameters  which  make  the  blast  properties  and  distances  dimen¬ 
sionless.  In  particular,  overpressure  and  dynamic  pressure  are  normalized 
with  respect  to  the  ambient  air  pressure  Pq,  and  distance  is  scaled  with 
respect  to  the  quantity  (E/pq)^/^,  where  E  is  the  energy  yield  of  the 
charge.  As  was  mentioned  earlier,  E  is  proportional  to  the  charge  weight 
W,  so  that  for  working  purposes  (E/p^)^*^^  can  be  replaced  by  (W/p^)^^^. 
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Table  1.  Ekjulvalent  we'.ght  factors  for  comparing  alt  blast  overpressures. 
Data  from  Swlsdak  (1975). 


Explosive 


Equivalent  weight 
relative  to  TNT 


Composition  A-3 
Composition  B 
Composition  C-4 
Cyclotoi  (RDX/TNT  70/30) 

HBX-1 

HBX-3 

H-6 

Minol  II 

Octol  (HMX/TNT  70/30  and  75/25) 
PETN 

Pentollte 

Pentollte  (high  overpressure) 

Picratol 

Tet  ryl 

Tetrytol  (Tetryl/TNT  75/25,  70/30, 

TNETB 

TNT 

Tritoi-al 

ANFO  (AN/FO  94/6) 

PBX-9404 

PBX-9010 

Nitroglycerin  Dynamite  (50%) 
Ammonia  Dynamite  (50%) 

Ammonia  Dynamite  (20%) 

Gelatin  Dynamite  (50%) 

Gelatin  Dynamite  (20%) 


1.09 

1.11 

1.37 
1.14 

1.17 

1.14 

1.38 
1.20 
1.06 

1.27 

1.42 

1.38 

0.90 

1.07 

65/35)  1.06 

1.36 
1.00 
1.07 
0.  82 

1.13 

1.29 

0.9 

0.9 

0.7 

0.8 

0.7 


Note:  Comparisons  made  for  overpressures  mostly  In  the  range  5-100 
Ibf/ln. ^  (0.03  to  0.7  MPa). 


With  Sachs  scaling,  characteristic  times  (e.g.  arrival  time,  phase 

duration)  are  scaled  with  respect  to  the  quantity  (E/pq)'-/^  (I/cq), 

where  Cq  is  the  acoustic  velocity  of  the  ambient  air.  Velocities  are 

?  1  /3 

scaled  with  respect  to  Cq.  Impulse  is  scaled  with  respect  to  (E/p^)  / 

cq.  In  all  cases,  distances  are  scaled  with  respect  to  (E/pq)^/^,  so 
that  scaled  distance  varies  with  the  ambient  pressure.  Sachs  scaling  ap¬ 
plies  over  the  range  of  distances  where  the  air  behaves  as  an  ideal  gas. 
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It  is  not  intended  for  use  very  close  to  the  explosion  source,  where  peak 

q  o 

overpressure  is  of  order  10  Ibf/in.  (7  MPa)  or  greater. 

Standard  values  of  pQ,  cq,  and  the  reference  values  of  density  and 
temperature  are  given  in  Figure  44. 

Reflection  and  Refraction  of  Airblast 

When  a  large  explosion  occurs  at  low  altitude,  the  spherical  shock 
front  travels  downward  and  outward  until  It  impinges  on  the  ground  surface 
and  reflects. 

The  reflection  of  acoustic  waves  and  weak  shocks  from  a  rigid  surface 
is  analogous  to  reflection  in  optics  (Fig*  6a).  Regular  reflection  occurs, 
with  the  angles  of  incidence  and  reflection  equal,  and  at  the  reflecting 
surface  the  pressure  of  the  incident  wave  doubles.  Natural  ground  surfaces 
are  often  porous  and/or  deformable,  permitting  some  wave  penetration  and 
energy  transfer  to  the  ground.  This  can  give  reflection  that  is  out  of 
phase  with  the  incident  wave. 

Directly  below  the  source  of  a  strong  shock  there  is  also  a  region  of 
regular  reflection,  but  the  angles  of  incidence  and  reflection  are  not 
equal,  and  the  reflected  wave  is  not  perfectly  spherical  (Fig.  6b).  The 


Incide.  ;  wove  Incident  Wove 

os  Reflection  Proqtesses  Hittmij  Surface 


a.  Regular  reflection  of  acoustic  waves  b.  Reflection  of  strong  shocks, 
and  weak  shocks.  with  Mach  stem  formation. 

Figure  6.  Reflection  of  shock  waves  at  a  rigid  surface. 


Figure  7,  Reflection  of  the  shock  wave  from  an  air  burst  and  formation  of 
the  Mach  stem. 


complications  arise  from  the  fact  that  the  air  behind  the  Incident  shock  is 
at  higher  pressure,  density  and  temperature  than  the  free  air  through  which 
the  incident  wave  traveled.  Thus  the  reflected  wave  travels  in  a  medium 
with  different,  and  spatially  varying,  properties.  As  a  result,  its  veloc¬ 
ity  and  pressure  are  affected.  The  pre-’sure  amplification  in  a  normally 
reflected  wave  is  much  more  than  the  factor  of  2  chat  applies  in  weak 
shocks  and  elastic  waves;  the  amplification  factor  could  be  in  the  range 
8  to  20  for  reflection  from  "hard"  surfaces  (Baker  1973;  see  also  Fig. 

15).  A  pressure  transducer  located  in  air  within  the  region  of  regular  re¬ 
flection  experiences  two  separate  shocks,  one  as  the  Incident  wave  passes 
"outward  bound,"  and  a  second  as  the  reflected  wave  returns. 

The  reflected  shock  travels  faster  than  the  Incident  shock  and,  in  the 
region  where  oblique  reflection  occurs,  the  reflected  shock  can  catch  up 
with  the  incident  shock,  fusing  with  it  to  form  a  single  front,  called  a 
Mach  front.  As  the  shock  front  expands,  the  size  of  the  fused  zone  in¬ 
creases  and  the  height  of  the  Mach  stem  Increases  with  distance  from  ground 
zero  (Fig.  7). 

The  locus  of  the  triple  point,  shown  schematically  in  Figure  7,  is 
given  quantitatively  for  a  range  of  charge  heights  in  Figure  8,  using 
scaled  variables. 


So  far,  the  air  has  been  assumed  homogeneous.  However,  for  large  ex¬ 
plosions  the  earth's  atmosphere  is  not  homogeneous,  but  rather  is  layered. 
The  effective  velocities  of  sound  waves  and  shock  waves  vary  with  altitude. 


V  '/ 

(m/kg  3)  (fr/ID5) 


ScalC'^  Hon^Ofttol  Rodius  to  Triple  Point 

Figure  8.  Height  of  triple  point  (or  Mach  stem)  as  a 
function  of  horizontal  radius  from  the  burst  axis  to 
the  triple  point.  Curves  are  given  for  various  charge 
heights.  All  linear  dimensions  are  scaled  with  respect 
to  the  cube  root  of  charge  weight. 


Multiplccation 


Figure  9.  Approximate  magni¬ 
tude  of  reinforcement  for 
ground  level  air  blast  fo¬ 
cused  at  long  range  by  atmo¬ 
spheric  conditions.  The 
graphs  represent  vertical 
profiles  of  wave  velocity. 

The  number  on  each  graph  is 
a  representative  molt Iplica- 
tlon  factor  for  the  type  and 
Intensity  of  the  velocity 
profile.  After  Baker  (1973). 


largely  because  of  changes  In  temperature  and  wind  speed.  From  eq  11, 
acoustic  velocity  Is  proportional  to  the  square  root  of  absolute  tempera¬ 
ture,  l.e.  c  =  (yR) T^/^.  Wind  does  not  affect  c  of  Itself,  but  It 
changes  the  velocity  relative  to  the  stationary  ground  surface  (see  eq  18 
below).  If  wave  velocity  Increases  with  altitude,  waves  are  refracted  and 
turned  back  towards  the  ground  surface.  This  can  be  thought  of  as  analo¬ 
gous  to  refr  .ctlon  and  internal  reflection  of  rays  In  optics.  The  rein¬ 
forcement  of  ground  level  air  blast  from  such  effects  Is  referred  to,  some 
times  loosely,  as  focusl.ig.  At  certain  "focus”  distances  from  an  explo¬ 
sion,  ground  level  air  blast  can  be  multiplied  by  factors  up  to  100  by  ver 
tlcal  wave  velocity  gradients  of  various  types  and  intensities  (Fig.  9). 
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The  atmospheric  factors  conttolllng  blast  can  best  be  appreciated 
by  considering  simplified  conditions.  In  the  unlikely  event  of  perfectly 
calm  conditions  with  no  vertical  gradient  of  temperature,  the  acoustic 
velocity  is  uniform  and  the  wave  front  from  an  explosion  propagates  in  a 
perfectly  spherical  pattern.  Thinking  in  terms  of  "rays"  that  are  ortho¬ 
gonal  to  the  wave  front,  this  condition  gives  straight  rays  radiating  out 
from  the  explosion  source  (Fig.  10a).  A  much  more  likely  situation  is  one 
where  the  temperature,  and  therefore  the  acoustic  velocity,  decreases  with 
increase  of  altitude  at  a  characteristic  gradient  (adiabatic  lapse  rate). 

In  this  case,  the  wave  front  is  distorted  from  the  perfect  spherical  pat¬ 
tern  and  the  "rays"  are  bent  in  an  upward  direction  (Fig.  10b).  This  con¬ 
dition  tends  to  funnel  shock  and  noise  upward,  reducing  overpressure  and 
noise  at  ground  level  for  long  ranges.  With  the  opposite  condition  (an  in¬ 
crease  of  temperature  and  acoustic  velocity  with  increasing  height)  the 
wave  front  speeds  up  as  it  rises,  causing  the  "rays"  to  bend  downward 
(Fig.  lOc).  This  increase  of  temperature  with  height,  known  as  a  tempera¬ 
ture  Inversion,  often  develops  in  the  lowest  layers  of  the  atmosphere  as  a 
result  of  overnight  cooling  in  calm  weather.  Air  blast  at  ground  level  can 
increase  considerably  when  a  temperature  Inversion  exists,  and  so  this  sit¬ 
uation  is  unfavorable  for  commercial  blasting  in  populated  areas.  Large 
explosions  give  significant  air  blast  effect  up  to  high  altitudes,  and  in 
these  circumstances  it  is  unlikely  that  there  will  be  a  uniform  temperature 
gradient  in  the  atmosphere.  With  discontinuities  and  reversals  in  the 
gradient  of  acoustic  velocity,  the  ray  paths,  interferences,  and  reflec¬ 
tions  can  become  complicated.  At  ground  level,  there  may  be  zones  of 
relative  silence,  and  zones  where  the  blast  is  intensified.  Figure  lOd 
gives  an  example  of  a  condition  where  the  velocity  gradient  reverses. 

Another  factor  is  wind,  which  has  both  direct  and  indirect  effects. 

The  direct  effect  is  that  sound  waves  travel  in  moving  air,  so  that  the  ef¬ 
fective  velocity  relative  to  the  ground  (c*)  is  the  vector  sum  of  acoustic 
velocity  c  and  the  wind  velocity  u: 

C'  =  C  +  U  CO80  (18) 

where  3  is  the  azimuth  from  the  downwind  direction,  u  varies  with  height, 
typically  because  of  wind  shear  against  the  ground  surface,  so  that  c'  also 
varies  with  height.  The  indirect  effect  comes  from  the  Influence  of  wind 
on  the  temperature  profile.  Direct  wind  effects  are  similar  to  temperature 
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a.  Undistorted  propagation  with  zero  gradient  of  acoustic 
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b.  Upward  deflection  when  acoustl*  velocity  decreases  with 
altitude. 
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Figure  10.  Refraction  of  air  waves  by  vertical  gradients  of 
acoustic  velocity. 
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d.  Deflection  and  Interference  with  stratification  and  reversal 
of  velocity  gradients. 

Figure  10  (cont'd). 


Wind  Oirection 


Figure  11.  Effect  of  wind  on  air  wave  propagation.  In  the  lower  atmosphere, 
rays  tend  to  be  bent  slightly  by  a  gradient  of  wind  velocity. 


effects,  given  corresponding  profiles  of  acoustic  velocity.  However,  while 
temperature  effects  are  omnidirectional,  wind  effects  vary  in  azimuth. 

With  a  typical  surfac  wind  profile,  rays  downwind  of  the  source  tend  to 
refract  downward,  while  rays  upwind  of  the  source  tend  to  bend  upwards. 
Figure  11  shows  how  blast  rays  would  be  affected  by  wind  alone,  i.e.  how 
they  would  leflect  with  wind  speed  increasing  with  altitude,  but  with  air 
temperature  uniform. 
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Lffect  of  Charge  Height,  or  Height  of  Burst 

VAien  a  charge  is  fired  close  to  the  surface,  che  attenuation  of  over¬ 
pressure  with  horizontal  distance  is  similar  to  the  attenuation  with  radius 
in  free  air.  However,  when  the  charge  1  high  above  the  surface,  the  over¬ 
pressure  experienced  at  ground  level  is  attenuated  by  both  downward  travel 
and  outward  travel,  i.e.  it  decreases  as  a  function  of  the  slant  range. 

This  naeans  that  as  the  horizontal  radius  from  ground  zero  decreases  to 
zero,  the  peak  overpressure  tends  to  a  limit  at  moderate  pressure  levels, 
Instead  of  increasing  to  very  high  values  as  it  would  if  the  charge  were 
close  to  the  surface.  As  the  charge  height  increases,  the  peak  ground 
level  overpressure  near  ground  zero  also  decreases.  However,  the  ground 
level  overpressure  at  large  scaled  radii  is  not  much  affected.  Figure  12 
gives  a  set  of  pressure  contours  which  shows  how  the  peak  overpressure  at 
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Figure  12.  Effect  of  charge  height  and  ground  level  radius  on 
peak  overpressure.  Contours  are  for  even  values  in  Ibf/ln.  ^ 

For  conversion  to  bars,  multiply  by  0.06895.  For  conversion  to 
kPa,  multiply  by  6.895.  Data  from  Swisdak  (1975). 
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Figure  13.  Effect  of  charge  height  and  ground  level  radius 
on  the  arrival  time  and  the  positive-phase  duration.  Data 
from  Swlsdak  (1975). 
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Figure  14.  Effect  of  charge  height  and  ground  level 
radius  on  the  positive  impulse.  The  curves  are  dravm 
for  given  levels  of  Impulse  in  (Ibf /in. ^)ms.  To  con¬ 
vert  to  (bars)tns  or  (kPa)ras,  multiply  the  values  by 
0.06895  or  6.895  respectively.  Data  from  Swisdak 
(1975). 


ground  level  varies  with  scaled  radius  from  ground  zero  and  with  scaled 
charge  height. 

The  arrival  time  and  the  positive  phase  duration  at  ground  level  are 
also  affected  by  the  charge  height.  The  general  pattern  of  variation  is 
shown  in  Figure  13. 

Corresponding  data  for  the  positive  impulse  are  given  iu  Figure  14. 


Attenuation  of  Air  Blast  and  Variation  of  Shock.  Front  Properties 

As  a  blast  wave  spreads  from  the  source,  It  attenuates.  The  peak 
overpressure,  the  dynamic  pressure,  the  shock  velocity  and  the  peak  parti¬ 
cle  velocity  all  decrease  with  increasl  'g  distance  from  the  source.  The 
density  and  temperature  of  the  air  behind  the  shock  front  also  decrease 
with  Increasing  radius.  The  positive  phase  duration  shows  an  overall  in¬ 
crease  with  Increasing  radius,  while  the  Impulse  of  positive  overpressure 
has  an  overall  decrease  as  radius  Increases.  Figures  15-18  Illustrate  the 
variation  of  shock  properties  with  radius  by  means  of  dimensionless  grapns 
based  on  Sachs  scaling. 


Figure  15.  Dimensionless  (Sachs) 
plot  of  peak  overpressure,  peak 
dynamic  pressure,  and  arrival 
time,  all  as  functions  of  radius. 
From  data  compiled,  selected  and 
adjusted  by  Baker  (1973). 
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Normolized  Density  and  Temperature  p/p  ,T/T^ 


Figure  16.  Dimensionless  (Sachs)  plot  of  shock  vel 
ocity  and  peak  particle  velocity  as  functions  of 
radius.  From  data  conpiled,  selected  and  adjusted 
by  Baker  (1973). 
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Figure  17.  Dimensionless  (Sachs)  plot  of  air  density  and  air 
temperature  as  functions  of  radius.  From  data  conplied,  se¬ 
lected  and  adjusted  by  Baker  (1973). 


Figure  18.  Dimensionless  (Sachs)  plot  of  positive 
phase  duration  and  positive  Impulse  as  functions 
of  radius.  From  data  compiled,  selected  and  ad¬ 
justed  by  Baker  (1973). 


Close  to  a  concentrated  source,  the  attenuation  of  overpressure  p^ 
should  be  rapid.  Taylor's  (1950)  solution  of  the  equations  for  spherical 
strong  shocks  predicts  an  inverse  cube  decay,  i.e.  p^  “  1/r^.  By  con¬ 
trast,  at  long  distances  from  the  explosion  the  air  blast  has  slowed  down 
to  become  an  acoustic  (elastic)  wave,  so  that  the  peak  pressure  should 
be  inversely  proportional  to  distance,  I.e.  p^  «  l/r.  Figure  19  gives  a 
representative  curve  for  the  variation  of  peak  overpressure  with  distance 
in  free  air.  To  avoid  complications  from  Inhoraogenelty  of  the  ambient  air 
and  from  surface  reflections,  it  is  necessarily  based  on  data  for  small 
charges  of  high  explosive,  which  have  finite  size.  The  broken  line  for 
high  pressure  at  close  range  represents  the  effect  of  finite  source  diam¬ 
eter;  for  a  true  point  source,  the  close-range  attenuation  should  follow  an 
inverse  cube  trend. 
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Figure  19.  Representative  attenuation  curve 
for  high  explosive  c'liatges  in  ''horuogeneous" 
free  air.  From  tabulated  unsealed  data  and 
calculated  extensions  of  curve,  as  given  by 
Swisdak  (1973).  Checked  against  other  field 
test  data  for  limited  ranges. 

In  typical  practical  problems  the  main  concern  is  with  air  blast  pro¬ 
pagated  at  ground  level  trom  a  charge  that  detonates  near  ground  level. 

The  overpressure  experienced  is  that  produced  after  reflection  and  wave 
fusion,  as  described  earlier.  To  remove  variability  due  to  charge  height 
(height  of  burst),  data  have  been  collected  for  surface  bursts  of  various 
sizes.  Test  results  are  not  completely  identical  when  presented  in  scaled 
form,  but  Figure  20  gives  a  representative  curve  for  the  attenuation  of 
peak  overpressure  with  scaled  distance  at  sea  level  (the  upper  limit  of  -he 
band  represents  a  1—kiloton  nuclear  burst  at  sea  level).  Comparing  Figures 
19  and  20,  it  can  be  seen  that  the  curves  for  free  air  and  for  surface  pro¬ 
pagation  are  not  much  different.  It  can  also  be  seen  that  the  data  do  not 
support  the  point-source  prediction  of  Inverse  cube  decay  (p^  «  l/r^). 

In  the  region  where  peak  overpressure  is  from  10  to  1000  Ibf/ln. ^  (70  kPa 
to  7  MPa),  an  inverse  square  decay  (p^j  «  1/r^)  is  a  better  approximation. 
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Figure  20.  Representative  attenuation  curve  for 
ground-level  overpressure  from  contact  and  near¬ 
surface  bursts.  HE  and  nuclear  data  from  various 
sources  (see  Swlsdak.  1975,  Baker  1973,  Glasstone 
1962,  Defense  Nuclear  Agency  1972,  U.S.  Army  1984, 
du  Pont  1977). 


Arrival  time  and  phase  duration.  Arrival  time  and  phase  duration  both 
scale  with  respect  to  the  cube  root  of  energy  or  charge  weight  (eq  15). 
Figure  21  gives  representative  curves  for  a  small  explosion  in  free  air. 
They  can  be  used  as  follows:  (1)  for  a  given  actual  distance,  obtain  the 
scaled  distance  by  dividing  by  the  cube  root  of  charge  weight,  (2)  for  this 
scaled  distance,  read  off  the  scaled  time,  (3)  convert  the  scaled  time  to 
actual  time  by  multiplying  by  the  cube  root  of  charge  weight. 


Figure  21.  Variation  of  arrival 
time  and  positive-phase  duration 
with  radius  for  an  explosion  In 
free  air.  From  Swlsdak  (1975). 


Air  blast  Impulse.  The  impulse  from  the  positive  phase  of  the  shock 
wave  Is  defined  by  eq  16  and  scaled  relations  for  Impulse  are  given  in  eq 
17.  As  a  very  rough  approximation,  the  pressure/t Ime  prot  for  the  positive 
phase  is  a  triangular  figure,  so  that  the  impulse  at  any  given  distance 
from  the  source  is  roughly  proportional  to  the  product  of  peak  overpressure 
(Fig.  20)  and  phase  duration  (Fig.  21).  Figure  22  gives  the  approximate 
magnitude  of  scaled  positive  iiiQ)ulse  as  a  function  of  scaled  radius  for 
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Figure  22.  Variation  of  positive  Impulse  with 
radius  for  surface  bursts  and  free  air  bursts. 
Adapted  from  Swlsdak  (1975). 
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surface  bursts  and  for  propagation  in  free  air.  Scaled  values  of  impulse 
are  converted  to  actual  values  by  multiplying  by  the  cube  root  of  charge 
weight. 

Dynamic  pressure.  If  the  air  behind  the  shock  front  behaves  as  an 
ideal  gas,  the  dynamic  pressure  q  can  be  expressed  as  a  function  of  the 
overpressure  p: 


q 


_ P _ 

(y-l)  +  2y  (po/p) 


(• 


_ P 

0.4  +  2. 8  (Pq/p) 


(19) 


where  Po  Is  ambient  pressure  and  Y  is  the  ratio  of  specific  heats,  usually 
taken  as  1.4  for  air.  If  this  relation  is  accepted,  the  attenuation  of 
peak  dynamic  pressure  q  is  determined  by  the  attenuation  relation  for  peak 
overpressure  p.  Equation  19  is  probably  acceptable  for  surface  bursts  when 
the  scaled  radius  exceeds  about  2  ft/lb^/^,  or  0.8  m/kg^/^. 

Particle  velocity.  The  particle  velocity  u  can  also  be  expressed  in 
terras  of  overpressure  p  and  ambient  pressure  Po  by  using  the  shock  equa¬ 
tions  and  the  gas  laws: 


u 


Co  P  Y+1  P  s-1/2 

Y  Po  ^  2y  Po^ 

0.714  cq  (p/pp) 

(.1  +  O.957(p/po)] 


(20) 


where  cq  Is  acoustic  velocity  in  ambient  air  ahead  of  the  shock  and  y  is 
taken  as  1.4. 

Shock  Velocity.  The  theoretical  relation  between  shock  velocity  U  and 
peak  overpressure  p  is 

U  .  c,  (H-  fi  • 

^  ^  2y  Pq-^ 

Cp  [l  +  0.857  (p/po)  (21) 

Actual  values  of  shock  velocity  can  be  obtained  from  test  data  for  arrival 
time,  as  shown  in  Figure  21.  The  shock  velocity  U  at  a  scaled  radius  r  is 

U  =  -^  (22) 

dt 


where  t  is  scaled  arrival  time  at  scaled  radius  r. 


Air  density  behind  the  shock  front.  Air  density  behind  the  shock, 
front,  p,  is  related  to  the  peak  overpressure  p  by 

P  _  1  +  I  (  Y+0/2y)p/po 

"pO  ”  1  +  [(y-1)/2y]p/Po 

1  +  0.857  p/po 

“  1  +  0.143  p/po  ^  ^ 

where  pq  is  air  density  ahead  of  the  shock. 

Reflected  pressure.  ITie  shock  equations  and  gas  laws  give  a  theoreti¬ 
cal  value  for  the  reflected  presure  p,.  at  the  Instant  when  the  incident 
shock  strikes  a  flat,  rigid  surface  at  normal  incidence: 

Pr  = 


where  p  is  the  Incident  peak  overpressure. 

Glasstone  (1962)  noted  that  for  very  high  values  of  the  incident  peak 
overpressure  p,  the  reflected  overpressure  tends  to  a  limit  of  Pr  =  8p» 
However,  Baker  (1973)  pointed  out  that  the  ideal  gas  assumption  is  probably 
not  justified  for  very  high  temperatures  and  pressures,  citing  studies  that 
suggest  Pf/p  can  be  as  high  as  20  or  more  for  very  Intense  shocks. 

Figure  23  shows  graphically  the  variations  of  dynamic  pressure,  parti¬ 
cle  velocity,  shock  velocity,  air  density  and  reflected  pressure  with  peak 
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Figure  23.  Variation  of  shock 
properties  with  peak  overpres¬ 
sure.  From  Swlsdak  (1975). 


overpressure.  The  range  of  p  Is  limited  In  the  graph  to  a  range  where  the 
Ideal  gas  assumption  Is  likely  to  be  justifiable. 

Air  Blast  from  Nuclear  Explosions 

In  general,  nuclear  explosions  are  on  a  bigger  scale  than  chemical  ex¬ 
plosions,  and  the  units  used  for  scaled  dimensions  are  based  on  a  kiloton 
of  yield*  rather  than  on  a  pound  or  kilogram  of  high  explosive.  The  very 
large  scale  of  the  blast  means  that  the  atmosphere  is  influenced  over  a 
large  range  of  altitude,  so  that  vertical  gradients  of  density  and  tempera¬ 
ture  are  important.  Another  difference  is  that  nuclear  explosions  emit  in¬ 
tense  thermal  radiation,  while  thermal  radiation  from  typical  high  explo¬ 
sives  is  negligible.  Thermal  radiation  absorbed  by  the  ground  surface 
causes  heating  and  the  formation  of  a  precursor  shock. 

A  nuclear  contact  burst  is  defined  as  a  nuclear  explosion  in  air  at 
scaled  heights  from  zero  to  5  ft/kt^/^  (1.5  m/kt^/^).  It  produces  a 
shock  front  that  Is  approximately  hemispherical.  An  err  burst  at  higher 
level,  say  from  5  to  160  ft/kt^/^  (1.5  to  50  m/kt^/-)  propagates  an  ap¬ 
proximately  spherical  shock,  which  reflects  on  contact  with  the  surface,  as 
described  previously  for  high  explosives  (p.  7).  A  Mach  stem  is  formed  by 
fusion  of  the  Incident  and  reflected  shocks,  and  air  motion  behind  the  Mach 
stem  produces  a  second  shock  (Fig.  24).  Ahead  of  the  Mach  stem,  heating  of 


MACH  STEM 


the  ground  by  thermal  radiation  causes  the  formation  of  a  precursor  wave 
near  the  surface  (Fig,  24).  The  amount  of  heating  produced  by  thermal  rad¬ 
iation  varies  with  the  reflectance,  or  albedo,  of  the  ground  surface,  and 
with  the  thermal  properties  of  the  surface  materials.  Ground  surfaces  are 
classified  as  either  "thermally  near-ideal"  or  "thermally  non-ideal";  the 
former  are  unlikely  to  create  a  precursor,  but  a  precursor  can  occur  over 
the  latter.  Table  2  gives  examples  of  these  two  types. 

Representative  distribution  of  peak  overpressure  is  shdwn  in  Figure  25 
for  surfaces  that  are  near-ideal.  Over  surfaces  that  are  thermally  non- 
ideal,  overpressure  contours  are  distorted  in  the  zone  of  Mach  reflection 
out  to  a  radius  of  1500  ft/kt^/^  (460  m/kt’^^),  or  to  the  6-lbf/ln. ^  con¬ 
tour. 


Scaled  Range  (m/kl''’) 
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a.  Very  high  overpressures  at  close  range. 

Figure  25.  Variation  of  peak  overpressure  with  charge  height  and  horizon¬ 
tal  radius  for  a  near-ideal  surface  subjected  to  a  nuclear  explosion  (U.S. 
Army  1984,  DNA  1972). 
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c.  Low  overpressures  at  appreciable  distances  from  the  explosion. 

Figure  25  (cont'd).  Variation  of  peak  overpressure  with  charge  height  and 
horizontal  radius  ^'or  a  near-ideal  surface  subjected  to  a  nuclear  explosion 
(U.S.  Array  1984,  DNA  1972). 
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d.  Very  Low  overpressures  in  the  far-field. 


Figure  25  (cont’d). 


Table  2.  Thermal  classification  of  ground  surfaces  (,U.S.  Army  1984), 


Thermally  near-ideal 
_ (precursor  unlikely) _ 

Water 

Ground  covered  by  white  smoke  layer 
Heat -ref lecting  concrete 
Frozen  tundra 
Ice 

Packed  snow 

[Moist]  Soil  with  sparse  vegeta¬ 
tion 

Commercial  and  industrial  areas 


Thermally  non-ideal 

_ (precursor  may  occur) _ 

Desert  sand 

Coral 

Asphalt 

Surfaces  wltn  low,  thick  vegetation 
Surfaces  covered  by  a  dark  smoke  layer 
Dark-colored  rock 

Most  agricultural  areas;  residential 
areas  in  cities 

Dry  soil  with  sparse  vegetation 


Attenuation  curves  for  air  blasts  of  nuclear  origin  are  given  in  Fig¬ 
ure  26.  One  curve  shows  the  variation  of  peak  overpressure  with  radius 
from  ground  zero  for  a  contact  burst.  The  other  curve  shows  the  variation 
of  peak  reflected  overpressure  with  slant  range.  It  gives  the  pressure  at 
ground  zero  as  a  function  of  burst  height,  or  the  pressure  as  a  function  of 
slant  range  for  points  within  a  30®  cone  under  the  burst  point. 
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Figure  26t  Attenuation  curves  for  air  blast 
from  nuclear  explosions.  The  solid  line  gives 
the  variation  of  peak  overpressure  with  ground- 
level  radius  for  a  contact  burst.  The  dashed 
line  gives  the  variation  of  peak  reflected  over¬ 
pressure  with  slant  range  for  an  air  burst  (for 
points  within  a  30“  cone  beneath  the  explosion). 

Data  from  U.S.  Army  (1964)  and  DNA  (1972). 

Dynamic  pressure  represents  a  wind  effect.  It  is  zero  at  ground  zero 
the  center  of  radial  symmetry;  it  increases  with  radius  to  a  maximum  and 
then  decreases  with  further  Increase  of  radius.  Dynamic  pressure  is  not 
affected  significantly  by  the  thermal  characteristics  of  the  ground  sur¬ 
face,  but  it  does  vary  with  the  amount  of  dust  that  is  raised. 
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Figure  27.  Variation  of  peak  dynaoilc  pressure  with  ground-level 
radius  for  a  nuclear  explosion  over  an  ideal  dust-free  surface 
(U.e  Army  1984,  DNA  1972). 


Figure  27  shows  the  variation  of  peak,  dynamic  pressure  with  radius  on 
an  ideal  surface  which  raises  no  dust.  At  close  range  the  dynamic  pressure 
is  Influenced  by  the  height  of  burst,  but  beyond  330  ft/kt^/^  (100 
m/kt^^^)  there  is  not  much  dependence  on  height  of  burst.  The  close- 
range  variation  of  peak  dynamic  pressure  with  charge  height  and  radius  is 
shown  in  Figure  28  for  an  ideal  surface. 

The  effect  of  dust  from  loose  surfaces  is  seen  in  the  zone  of  Mach  re¬ 
flection  (Fig.  29,  30).  The  effect  is  strongest  for  burst  heights  lower 
than  400  ft/kt^^^  (120  m/kt'’>^^),  and  at  ranges  closer  than  1200  ft/kt^/^ 
(370  m/kt^'^^).  Between  400  and  1200  ft/kt^^^  from  ground  zero,  dust  in¬ 
creases  the  peak  dynamic  pressure  in  comparison  with  dust-free  conditions 
(Fig.  27,  28).  At  closer  range,  dust  can  lower  the  peak  dynamic  pressure. 

The  Impulse  for  the  positive  phase  is  the  integral  of  overpressure 
with  respect  to  time,  and  so  it  is  affected  by  the  thermal  properties  of 
the  ground  surface.  Figure  31  gives  a  representative  distribution  of  im- 
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Figure  28.  Variation 
of  peak,  dynamic  pres¬ 
sure  with  charge  height 
and  radius  for  a  nu¬ 
clear  explosion  over 
an  ideal  dust-free  sur¬ 
face  (U.S.  Army  1984, 
DNA  1972). 
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Figure  29.  Variation  of  peak,  dynamic  pressure  with  charge  height  and 
radius  for  a  nuclear  explosion  over  a  surface  where  light  dust  condi¬ 
tions  are  generated  (U.S.  Army  1984,  DNA  1972). 
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Figure  30.  Variation  of  peak  dynamic  pressure  with  radius  for  nu¬ 
clear  air  bursts  over  surfaces  which  generate  heavy  dust  conditions. 
Adapted  from  U.S.  Array  (1984)  and  DNA  (1972). 
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Figure  31.  Variation  of  positive  impulse  with  charge  height  and  radius 
for  nuclear  explosions  over  a  near-ideal  surface  (U.S.  Army  1984  DNA 
1972). 
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FigucG  3''!.  Variation  of  positive  phase  duration  (overnressure)  with  charge 
height  and  radius  for  nuclear  explosions  over  near-ideal  surfaces  (U.S. 

Army  1984,  DNA  1972). 

pulse  for  a  surface  that  Is  near-ideal,  and  Figure  32  shows  how  the  posi¬ 
tive  impulse  at  the  surface  attenuates.  The  distribution  changes  sotnewhat 
when  the  surface  Is  non-ideal.  Figure  33  shows  how  the  attenuation  curves 
for  positive  impulse  at  the  surface  are  affected;  in  general,  the  Impulse 
at  close  range  tends  to  Increase  slightly  by  comparison  with  the  near-ideal 
surface. 

Ttie  distribution  of  phase  duration  for  positive  overpressures  is  shown 
in  figure  34  for  a  surface  that  has  near-ideal  thermal  properties.  When 
the  surface  is  non-ideal,  the  distribution  is  distorted  significantly,  even 
for  very  high  bursts  and  for  large  distances  from  ground  zero  (see  U.S. 

Army  Technical  Manual  5-858-2  for  details). 

The  phase  duration  for  positive  dynamic  pressure  is  shown  in  Figure 
35.  This  distribution  is  for  surfaces  where  little  or  no  dust  is  kicked  up 
by  the  blast.  When  heavy  dust  is  raised,  phase  duration  changes  signi¬ 
ficantly  (see  U.S.  Army  Technical  Manual  5-858-2  for  details). 
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Figure  35.  Variation  of  positive  phase  duration  (dynamic  pressure) 
with  charge  height  and  radius  for  nuclear  explosions  over  surfaces 
where  little  or  no  dust  Is  generated  (U.S.  Army  1984,  DNA  1972). 


The  decay  of  pressure  waves  with  time  at  a  fixed  location  has  been  de¬ 
scribed  mathematically  by  variants  of  the  function  xe""^.  Both  overpres¬ 
sure  and  dynamic  pressure  have  their  maximum  values  at  the  shock  front  and, 
for  a  given  location,  both  decrease  to  zero  at  about  the  same  time  (the 
positive  phase  of  the  dynamic  pressure  is  actually  slightly  longer  than 
that  of  the  overpressure).  Except  in  regions  of  very  high  overpressure 
Cpm  ^  Ibf/ln.  or  >  0.  5  MPa),  peak  dynamic  pressure  is  less  than  peak 
overpressure,  and  their  rates  of  decay  are  therefore  different.  The 
exponential  decay  of  overpressure  and  dynamic  pressure  at  a  given  point  has 
been  described  by  the  empirical  relations 

^  exp  ,  '  (25) 

^  =  (1  -  exp  (-  (26) 

where  p  and  q  are  the  tra  -lent  overpressure  and  dynamic  pressure  at  time 
t,  and  t+  is  the  positive  phase  duration  (Glasstone  1962).  Figures  36  and 
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Normalized  Time 


Curves  from  U.S.  Army  (1984)  and  DNA  (1972). 


b.  Curves  from  Glasstone  (1962). 


Figure  37.  Idealized  exponential  decay  of  positive 
dynamic  pressure  with  time  at  sea  level.  See  text 
for  normalizing  parameters. 


37  give  some  accepted  approximations  of  the  wave  decay  for  overpressure  and 
dynamic  pressure. 


Air  Blast  from  Underground  Explosions 

Air  blast  is  muffled  when  the  charge  is  detonated  below  the  ground 
surface.  In  general,  there  are  two  identifiable  pulses  in  the  air  blast 
record  from  an  underground  explosion.  The  first  is  identified  with  ground 
shock,  i.e.  there  is  an  Initial  transfer  of  energy  to  the  air  when  the 
shock  wave  hits  the  ground/air  interface.  The  second  is  induced  by  ground 
motion  and  venting  as  the  gas  bubble  bursts  through  the  surface.  In  hard 
ground  that  can  transmit  shock  without  much  dissipation,  the  first  pulse. 
Induced  by  ground  shock,  may  be  dominant.  In  soft  ground  that  absorbs 
shock,  the  second  pulse,  induced  by  venting,  may  be  the  dominant  one. 

Estimates  of  air  blast  from  burled  charges  are  based  on  the  data  for 
surface  charges  (Fig.  20,  26).  An  air  blast  prediction  is  made  for  a  sur¬ 
face  explosion,  with  appropriate  adjustments  for  ambient  environmental  con¬ 
ditions.  Thereafter,  the  predicted  overpressure  can  be  obtained  either  by 
applying  a  danping  factor  to  the  surface-burst  pressure,  or  by  applying  a 
multiplying  factor  to  the  range  scale  of  the  curve  for  surface  bursts. 

In  the  method  given  by  Johnson  (1971),  the  predicted  peak  overpressure 
for  a  surface  burst  is  multiplied  by  a  damping  factor  (Fig.  38),  and  the 
result  is  the  peak  overpressure  for  the  dominant  pulse  of  the  underground 
explosion,  without  regard  for  whether  the  ground  shock  pulse  or  the  venting 
pulse  is  the  dominant  one.  The  dandling  factor  is  a  function  only  of  charge 

Scaled  Chorge  Depth  ((t/lb"^) 


Figure  38.  Air  blast  damp¬ 
ing  factor  as  a  function  of 
charge  depth  for  large  bur¬ 
ied  explosions.  The  air 
blast  overpressure  for  a  sur¬ 
face  charge  is  multiplied 
by  the  appropriate  damping 
factor  to  predict  the  air 
blast  overpressure  for  a 
buried  explosion.  After 
Johnson’s  (1971)  adaptation 
of  a  relation  by  Reed. 
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Figure  39.  Range  multiplication  factor  K  for 
predicting  air  blast  from  burled  explosions. 

To  obtain  the  overpressure  p  at  a  specified 
scaled  range  F  from  the  underground  explosion, 
read  from  Figure  16  (or  equivalent  curve)  the 
value  of  p  at  scaled  radius  Kr.  Graphical 
version  of  computational  procedure  given  by 
Swlsdak  (1975). 


depth,  without  reference  to  ground  material.  In  the  method  described  by 
Swlsdak  (1975),  the  range  is  adjusted  by  a  multiplying  factor,  such  that 
the  range  for  a  given  overpressure  becomes  greater  for  a  burled  explosion 
than  for  a  surface  explosion.  The  multiplication  factor,  which  is  a  func¬ 
tion  of  charge  depth  and  ground  density,  is  obtained  by  raising  e  to  a  di¬ 
mensional  power.  This  technical  impropriety  could  be  remedied  by  convert¬ 
ing  scaled  charge  depth  in  ft/lb*-/^  to  a  dimensionless  value  from  Figure 
5,  but  for  present  purposes  it  is  easier  to  give  the  multiplier  K  in  graph 
icai  form  (Fig.  39).  To  predict  the  overpressure  p  at  a  specified  range  r 
for  a  given  underground  explosion  (weight  W),  an  appropriate  value  of  K  is 
obtained  from  Figure  39,  r  is  scaled  as  r/W^/3,  and  p  is  read  from  a  sur¬ 
face  burst  curve  (Fig,  20)  for  a  distance  of  K  r/W^/^, 

When  a  row  of  cratering  charges  is  fired  simultaneously  the  peak  over 
pressure  at  a  distance  is  usually  reinforced.  Figure  40  gives  multlpllca- 


Figure  40.  Overpressure  multiplier 
for  air  blast  from  a  rcw  charge.  If 
the  point  of  concern  is  in  line  with 
row  of  charges,  use  the  lower  line. 
If  the  point  of  concern  is  off  to 
one  side  of  the  row,  and  reasonably 
distant,  use  the  upper  line.  From 
Johnson  (1971). 


tlon  factors  for  row  charges.  This  is  based  on  observations  for  row 
charges  at  practical  spacings.  In  principle,  there  are  charge  spacings  and 
range  positions  which  do  not  necessarily  give  amplification;  waves  from  the 
Individual  charges  can  interfere  in  a  variety  of  ways,  including  superposi¬ 
tion  of  positive  and  negative  phases. 

In  normal  blasting  operations,  air  blast  is  not  the  controlling  damage 
criterion;  ground  motion  is  likely  to  become  a  problem  before  air  blast 
causes  damage. 

Air  Blast  from  Underwater  Explosions 

Detonation  underwater  greatly  reduces  airblast  in  con^jarison  with  sur¬ 
face  detonation,  even  at  very  shallow  scaled  depth.  Figure  4 1  gives  con¬ 
tours  of  peak  overpressure,  showing  how  air  blast  varies  vdth  charge  depth 
and  with  horizontal  distance  from  "surface  zero,"  the  point  directly  above 
Che  charge.  Contours  are  given  for  horizontal  planes  at  three  different 
heights  above  the  water  surface.  Point  values  of  peak  overpressure  for  a 
surface  burst  over  solid  ground  are  also  shown  on  the  graph*.  Data  for 
greater  scaled  heights  are  given  by  Swlsdak  (1978). 

*  Attenuation  curves  given  by  Hallanger  (1976)  for  air  blast  overpressure 
from  underwater  explosions  seem  questionable.  The  indicated  air  blast  from 
some  shallow  underwater  explosions  is  greater  than  air  blast  from  surface 
bursts,  and  the  indicated  attenuation  is  acoustic  (n  =  -1)  for  all  levels 
of  overpressure. 
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Peak  overpressure 
for  surface  burst  over  solid  ground 
at  sea  level  (Ibf/in^) 


Figure  41.  Mr  blast  overpressure  from  under¬ 
water  explosions.  The  plot  shows  the  effect  of 
charge  depth,  radius  from  surface  zero,  and  the 
height  of  observation  above  water  level.  Fig¬ 
ure  Ala  shows  the  general  pattern,  while  Figure 
Alb  gives  more  detail  for  fairly  long  ranges. 
Intersections  of  the  contours  with  the  axis  for 
zero  charge  depth  give  the  attenuation  for  a 
surface  burst  over  water.  The  figures  in  par¬ 
entheses  are  the  surface-level  overpressure 
values  that  occur  at  the  indicated  distances 
for  a  surface  burst  over  solid  ground.  Adapt¬ 
ed  from  Swisdak  (1975). 


A2 


Structural  damage  from  air  blast  depends  on  the  characteristics  cf  the 
air  shock,  principally  peak  overpressure  and  phase  duration,  and  also  on 
the  characteristics  of  the  structural  element,  mainly  frequency  response 
and  strength.  In  cases  of  severe  damage  where  dynamic  pressure  is  a  fac¬ 
tor,  the  drag  characteristics  and  "porosity"  of  the  target  may  also  be  im¬ 
portant. 

Damage  criteria  are  based  on  peak  overpressure,  making  the  tacit  as¬ 
sumption  that  phase  duration  does  not  vary  significantly  for  a  given  over¬ 
pressure.  Such  an  assumption  seems  justifiable  for  limited  ranges  of 
explosive  yield,  since  phase  duration  follows  cube  root  scaling,  but  signi¬ 
ficantly  different  waveforms  are  possible  (e.g.  the  long-period  N-wave  of  a 
sonic  boom).  Approximate  overpressure  levels  for  positive  damage  are  shown 
in  Figure  42,  together  1th  the  safe  limit  accepted  by  the  U.S.  Bureau  of 
Mines. 
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Figure  42.  Approximate  values  of  overpres¬ 
sure  at  which  various  kinds  of  damage  are 
likely  to  occur.  Data  from  various  sources, 
e.g.  Glasstone  (1962),  Siskind  et  al.  (1980), 
du  Pont  (1977). 


In  civil  blasting  operations  the  emphasis  is  on  limiting  air  blast  to 
safe  or  tolerable  levels,  and  on  verifying  the  limits.  Air  blast  at  these 
levels  is  commonly  perceived  as  noise  of  relatively  high  frequency  rather 
than  as  shock,  and  safe  limits  are  usually  expressed  In  decibels  (dB). 
Measured  values  depend  to  some  extent  on  the  response  characteristics  of 
the  sensor,  and  safe  limits  are  specified  with  reference  to  the  measuring 
device.  Table  3  gives  the  USBM  safe  limit  for  various  types  of  measuring 
systems. 

Table  3.  USBM  safe  limits  for  air  blast. 

Safe  limit 


-  m-  II  O 

Characteristic  of  measuring  system  dB  Ibf /in. 


0. 1  Hz  high-pass  134  0.0145 

2  Hz  high-pass  133  0.013 

5  Hz  or  6  Hz  high-pass  129  0.0082 

C-slow  (events  not  exceeding  2  s  105  0.00052 

duration) _ 


Air  blast  intensity  in  decibels  I  is  related  to  the  pressure  level  p 
by 

I  =  20  logio  (p/p*)  dB  (27) 

where  p*  is  a  reference  level  of  air  pressure  taken  as  2x10"**  dyne/cm^ 

(2. 9x10"^  Ibf/in. ^)  at  20®C. 

Effects  of  Ambient  Pressure  and  Temperature 

The  effects  of  an  explosion  in  free  air  are  much  the  same  in  cold  re¬ 
gions  as  they  are  anywhere  else,  although  the  temperature  and  pressure  of 
the  ambient  air  can  produce  some  minor  quantitative  differences.  Cold  re¬ 
gions  tend  to  have  lower  air  temperatures  and  some  regions,  notably  Green¬ 
land  and  Antartlca,  have  high  surface  elevations,  with  consequent  low 
values  of  atmospheric  pressure  and  density.  Figures  43  and  44  show  how 
pressure,  temperature,  density  and  acoustic  velocity  vary  with  altitude  in 
free  air. 

In  a  homogeneous  air  mass,  the  peak  air  blast  overpressure  p^  at  a 
given  scaled  radius  varies  with  the  ambient  air  pressure  pg  such  that  the 
shock  strength  is  constant  for  constant  scaled  radius.  Since  the  shock 
strength  is  the  ratio  of  absolute  pressure  behind  the  shock  front  to  the 
ambient  pressure,  t.e.  (Pm+Po)/?  o>  write: 
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+  1 


(28) 
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P  01 


!s2 

P02 


where  subscripts  1  and  2  denote  two  different  levels  of  ambient  pressure 
(or  altitude),  and  the  scaled  radius  is  constant.  Thus,  at  a  given  scaled 
distance  from  an  explosion,  the  peak  overpressure  will  vary  directly  with 
the  atmospheric  pressure  as 


^in2  ^  p^ 

Pml  "  POl 


(29) 


In  other  words,  if  p^i  is  the  peak  overpressure  given  by  standard  curves 
adjusted  to  sea  level  pressure  (I  atmosphere,  or  14.7  Ibf/in. ^),  that  fig¬ 
ure  has  to  be  adjusted  by  a  factor  (po?./Poi)  obtain  the  peak  overpres¬ 
sure  at  the  same  scaled  radius  when  the  atmospheric  pressure  is  Po2* 
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Figure  43.  Variation  of  atmospheric  pro¬ 
perties  with  altitude  according  to  U. S. 
Standard  Atmosphere  of  1962. 
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Figure  44.  Relative  values  of  atmo 
spheric  properties  as  functions  of 
altitude,  according  to  U.S.  Stand¬ 
ard  Atmosphere  of  1962. 


The  effect  of  atmospheric  pressure  on  the  scaled  radius  r  for  a  fixed 
value  of  Pm/ro  giveu  by  Sachs  scaling  as! 


-  "  ^01^ 


where  rj  and  t2  are  the  scaled  radii  where  a  fixed  value  of  Pm/Po  Is 
experienced  with  ambient  pressure  levels  pgj  and  Po2‘ 

The  net  effect  produced  by  a  change  of  atmospheric  pressure  is  esti¬ 
mated  from  a  two-step  procedure.  The  scaled  radius  is  adjusted  in  accord¬ 
ance  with  eq  30,  and  the  peak,  overpressure  for  that  adjusted  scaled  radius 
at  the  reference  elevation  is  multiplied  by  the  factor  given  by  eq  29. 
Figures  15-18  indicate  how  £imblent  pressure  influences  the  standard  blast 
parameters. 

Ambient  air  temperature  affects  density,  wave  velocity,  arrival  time 
and  phase  duration.  Combining  the  effects  of  atmospheric  pressure  and  am¬ 
bient  air  temperature,  the  initial  density  varies  as 


Sonic  velocity  In  the  ambient  air  varies  as 


^02 


1/2 


(32) 


and  shock,  velocity  varies  in  th<  same  way. 
adjusted  radius  r2  given  by  eq  JO  is 


^  (I  1  1/3  Ti  1/2 

ti  '■p  2^ 


The  arrival  time  t2  at  the 


(33) 


and  positive  phase  duration  varies  in  the  same  way.  The  Impulse,  which  is 
the  integral  of  overpressur i  with  respect  to  time,  is  adjusted  by  a  factor 
made  up  from  the  adjustment  factors  for  overpressure  and  time,  i.e. 


.PQ2..2/3  ,T^,.i/2 
'^POI^  ^T2-’ 


(34) 


Convenient  tables  and  graphs  for  adjusting  blase  parameters  in  accord¬ 
ance  with  Sachs  scaling  are  given  by  Swisdak  (1975). 

Explosions  in  Vacuum  or  in  Space 

Having  noted  that  the  standard  altitude  adjustments  are  not  valid  for 
very  big  differences  of  aicitude*  it  is  wort’u  considering  briefly  the  lim¬ 
iting  condition  when  pressure  tends  to  zero.  In  a  vacuum,  or  in  space,  the 
ambient  pressure  and  density  are  effectively  zero.  If  the  explosion  source 
functions  normally,  it  generates  explosion  products  at  high  initial  pres¬ 
sure  and  temperature,  but  it  cannot  impart  mechanical  disturbance  to  the 
surroundings  because  there  is  nothing  there.  Since  the  explosion  products, 
in  the  form  of  particles  and  gas  molecules,  receive  momentum  from  the  ex¬ 
plosion  process,  they  ought  to  travel  outward  indefinitely  without  resist¬ 
ance,  but  subject  to  gravity.  If  the  explosion  source  produces  radiation, 
as  a  nuclear  source  does,  the  radiation  will  be  transmitted,  subject  only 
to  geometrical  attenuation. 

Air  Blast  Attenuation  Over  Snow  Surfaces 

Air  blast  measurements  have  been  made  for  propagation  over  the  surface 
of  deep  snow,  and  some  results  are  given  in  Figure  45.  Since  all  these 
measurements  were  made  at  high  altitude,  the  data  are  adjusted  to  a  common 
surface  elevation  of  6000  ft  (1830  m).  For  conparison  with  standard  sea 
level  data,  the  curve  of  Figure  45  has  to  be  shifted  up  in  accordance  with 
the  adjustment  procedures  described  earlier.  At  overpressures  above  1 
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Figure  45.  Air  blast  attenuation  over  the  surface  of  very  deep 
snow-  All  measurements  were  made  at  high  altitude  and  the  data 
have  been  adjusted  to  a  common  surface  elevation  of  6000  ft 
(1830  m). 
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Ibf/in.  (7  kPa)  there  Is  close  agreement,  but  the  Ear-field  measurements 
over  snow  are  lower  than  the  values  predicted  by  extrapolation  of  the 
standard  curves.  The  rate  of  attenuation  for  the  Ear  field  is  such  that 
overpressure  is  approximately  proportional  to  r"^*^. 

Measurements  of  overpressure  and  positive-phase  impulse  were  also  made 
over  seasonal  snow  cover  (Wlsotskt  and  Snyer  1966),  using  small  charges  (1 
lb,  or  0.45  kg)  and  scaled  thicknesses  for  the  snow  cover  in  the  range  0.23 
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Figure  46.  Air  blast  attenuation  over 
seasonal  snow  (scaled  thickness  0.23 
to  1.65  m/kg^/3^  or  0.58  to  4.2  ft/ 
lbl/3),  Th.3e  measurements  were  made 
at  an  altitude  of  10,800  ft  (3292  m) . 
Summary  of  results  by  Wisotski  and 
Snyer  (1966). 
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Figure  47.  Attenuation  of  positive  impulse  over  sea¬ 
sonal  snow  (thickness  and  site  altitude  as  for  Fig. 
46).  Summary  of  results  by  Wisotski  and  Snyer  (1966). 


to  1.65  m/kgl/3  (0.53  to  4.2  ft/lb^'^^).  The  altitude  of  the  site  was 
3292  m  (10,800  ft),  with  average  ambient  pressure  0.68  bar.  The  oversnow 
attenuation  curves  for  overpressure,  with  charge  height  as  parameter,  are 
shown  in  Figure  46.  For  the  range  of  the  data,  overpressure  is  approxi- 
mately  proportional  to  1/r  . 

The  corresponding  Impulse  data  of  Wisotski  and  Snyer  show  considerable 
scatter;  the  trend  lines  are  given  in  Figure  47.  Impulse  is  roughly  pro¬ 
portional  to  ^  for  this  range. 

Shock  Reflection  from  Snow  Surfaces 

Simple  theory  for  shock  reflection  assumes  a  rigid  reflecting  surface, 
so  that  no  energy  is  lost  in  the  reflection  process.  However,  when  the  re¬ 
flecting  surface  can  absorb  energy  by  inelastic  deformation,  the  amplitude 


Figure  48.  Reflected  pressure  as  a 
function  of  Incident  pressure  for  rel¬ 
atively  strong  air  blast  over  deep  ice 
cap  snow  (measured)  and  over  a  rigid 
surface  (theoretical).  After  Ingram 
(1962), 


Figure  49.  Height  of  the  triple  point  (Mach  stem)  as 
a  function  of  radius  and  charge  Height  for  shock  re¬ 
flection  over  deep,  dense  snow  at  high  altitude.  From 
Ingram  (1962). 


of  the  reflected  wave  ought  to  be  less  than  that  for  reflection  from  a  rig¬ 
id  surface. 

Figure  48  summarizes  the  results  of  measurements  made  on  the  Greenland 
ice  cap,  where  air  blast  from  high  explosive  charges  was  reflected  from  the 
surface  of  dense,  but  very  deep,  snow.  The  reflected  pressure  was  signifi¬ 
cantly  lower  than  the  theoretical  reflected  pressure  over  a  rigid  surface, 
as  calculated  from  eq  23. 

Figure  49  gives  some  data  for  the  height  of  the  triple  point  (heiglit 
of  Mach  stem)  for  shock  reflection  over  deep  snow  at  high  altitude.  Com- 
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paring  these  curves  with  the  "standard"  curves  of  Figure  8,  It  can  be  seen 
that  the  snow  surface  curves  are  displaced  downward.  In  other  words,  the 
locus  of  the  triple  point  is  closer  to  the  surface  than  might  be  expected. 

Hisotski  and  Snyer  (1966)  made  measurements  over  seasonal  snow  cover 
(thickness  0.23  to  1.65  m/kg^/^)  and  also  over  bare  ground.  Their  mea¬ 
surements  of  peak  overpressure  p^  for  the  region  of  Mach  reflection  show 
that  Pm  is  higher  over  bare  ground  than  over  snow  for  given  values  of 
radius  and  charge  height  (Fig.  50).  The  general  pattern  of  behavior  is 
also  illustrated  in  Figure  51. 


Scoi»d  Rodiui  from  Ground  Z.ro  Scoldd  Rodiuo  from  Ground  Ztro 

Figure  50.  Peak  reflected  overpressure  (Mach  region)  as  a  function  of  radius 
and  charge  height,  conparlng  measurements  over  seasonal  snow  and  measurements 
over  bare  ground  (see  Fig.  46  for  site  conditions).  Summary  of  results  by 
Wisotskl  and  Snyer  (1966). 
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Figure  52.  Height  of  triple  point  (Mach  stem)  as  a  function  of  radius  and 
charge  height  for  explosions  over  seasonal  snow  cover  and  over  bare  ground 
at  high  altitude.  From  Wisotski  and  Snyer  (1966). 

The  path  of  the  triple  point  over  seasonal  snow  and  over  bare  ground 
is  shown  in  Figures  50  and  52.  In  general,  the  Mach  stem  is  higher  over 
bare  ground  than  over  snow  at  the  same  radius. 

Wisotski  and  Snyer  defined  a  reflection  coefficient  K  as  the  ratio  of 
charge  weights  which  will  produce  equal  overpressures  at  a  given  radial 
distance  when;  (a)  th.e  charge  is  fired  in  free  air,  and  (b)  the  charge  is 
fired  close  to  the  Mach  reflection  surface.  Designating  these  charge 
weights  as  and 

K  =  Wfa''Wmr  (35) 

If  W|^{^  is  taken  as  unit  weight,  K  can  be  expressed  in  terras  of  scaled 
radii  to  a  given  contour  of  peak  overpressure: 


(36) 


where  rp^  is  scaled  radius  in  free  air  to  peak  incident  overpressure 
level  p^j,  and  r^R  is  scaled  radius  to  peak  reflected  pressure  p^  in 
the  Mach  region.  Values  of  K  for  bare  ground,  snow  and  concrete  are  plot¬ 
ted  against  scaled  charge  height  in  Figure  53. 
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Figure  53.  Reflection  coefficient  (see  text)  as  a  function 
of  charge  height  for  snow,  bare  ground  and  concrete  at  high 
altitude.  From  Wlsocskl  and  Snyer  (1966). 
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Figure  54.  Positive  impulse  as  a  function 
of  radius  and  charge  height  for  explosions 
over  seasonal  snow  cover  and  over  bare 
ground  at  high  altitude.  Summary  of  re¬ 
sults  by  Wlsotskl  and  Snyer  (1966). 
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Figure  54  (cont'd).  Positive  impulse  as  a 
function  of  radius  and  charge  height  for  ex¬ 
plosions  over  seasonal  snow  cover  and  over 
bare  ground  at  high  altitude.  .  Summary  of 
results  by  Wisotskl  and  Snyer  (1966). 
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Figure  54  (cont'd). 


Impulse  measurements  over  seasonal  snow  and  over  bare  ground  gave 
fairly  scattered  data;  trend  lines  ate  shown  In  Figure  54.  In  general,  the 
Impulse  over  bare  ground  Is  higher  than  the  Impulse  over  snow. 

Shock  reflection  from  a  surface  of  land-based  Ice  should  not  be  signi¬ 
ficantly  different  from  rigid-surface  reflection.  However,  reflection  from 
a  floating  Ice  cover  could  Involve  energy  dissipation  and  reduction  of  amp¬ 
litude  In  the  reflected  wave.  At  t.ils  time,  no  test  data  are  available. 

Shock  Velocity  Over  Snow 

Iteasurements  of  shock  velocity  were  made  on  the  Greenland  ice  cap 
(Fig.  55; .  They  agree  fairly  well  with  the  theoretical  shock  velocity,  as 
calculated  from  eq  20.  The  discrepancies  can  probably  be  attributed  to 
calibration  errors  in  the  pressure  transducers,  and  perhaps  to  other 
problems  in  the  measuring  equipment  (one  set  of  data  shows  shock  velocity 
dropping  below  the  expected  acoustic  velocity). 
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Figure  55.  Shock  velocity  as  a  function  of  peak  over¬ 
pressure  for  explosions  over  deep  snow  at  high  alti¬ 
tude.  After  Smith  (1962). 


Variation  of  Shock  Pressure  with  Charge  Height  Over  Snow 

At  ground  zero,  the  point  on  the  surface  directly  below  an  air  burst, 
peak  overpressure  decreases  as  charge  height  increases,  in  accordance  with 
the  free-air  attenuation  relation.  For  a  given  charge  height  R,  overpres¬ 
sure  decreases  with  horizontal  radius  from  ground  zero  (R)  in  accordance 
with  the  increase  of  slant  range  S  from  the  charge,  where  S  = 

Figure  56  gives  contours  of  peak  overpressure  for  air  blast  over  deep  snow, 
In  the  same  form  that  Is  used  for  standard  curves  (see  Fig.  12,  25).  Fig¬ 
ure  57  gives  the  data  In  another  form,  showing  attenuation  curves  with 
charge  height  as  parameter. 

The  results  of  Wisotski  and  Snyer  (1966)  give  attenuation  curves  mea¬ 
sured  over  seasonal  snow.  Curves  are  given  for  a  range  of  charge  heights 
(Fig,  46),  and  a  free-air  attenuation  curve  Is  Included  for  comparison.  At 
any  given  radius  from  ground  zero,  the  Mach  region  peak  overpressure,  p^, 
decreases  as  charge  height  increases,  within  the  limited  range  of  the  data. 
The  results  are  shown  in  different  form  in  Figure  51,  which  does  not  dis¬ 
play  the  characteristic  "bulges"  seen  in  Figures  12  and  56  for  somewhat 
broader  fields  of  coverage.  For  any  finite  charge  height,  Pj^  is  higher 
than  the  free-air  value  after  Mach  reflection  from  snow.  However,  for  zero 
charge  height,  with  the  charge  lying  in  the  snow,  p^j  is  smaller  than  the 
free-air  value. 
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Release  of  Avalanches  by  Air  Blast 

Snow  lying  on  a  slope  can  slip  and  form  an  avalanche  when  the  down- 
slope  shear  stress  exceeds  the  shear  strength  at  any  level  In  the  snow 
pack.  Under  natural  circumstances,  such  Instability  develops  when  stress 
Increases  faster  than  strength  because  of  rapid  accumulation,  or  alterna¬ 
tively  when  snow  strength  deteriorates  under  constant  stress.  If  a  snow 
slope  Is  highly  stable  (gentle  slope,  well-settled  dry  snow)  It  is  unlikely 

to  be  disturbed  by  air  blast.  However,  If  the  snow  Is  only  marginally 

stable  (slope  angle  30°-60'’,  stress/strength  ratio  approaching  unity),  then 
air  blast  could  trigger  an  avalanche.  The  question  Is,  how  much  overpres¬ 
sure  Is  required,  and  over  how  wide  an  area  must  it  be  distributed? 

Avalanche  folklore  includes  instances  of  avalanche  release  by  the 
sound  of  church  bells,  alpenhorns  or  bellowing  voices,  but  it  is  hard  to 

believe  that  such  weak  elastic  waves  have  much  effect.  At  the  other  ex¬ 

treme,  marginally  stable  snow  Is  very  likely  to  be  released  by  detonation 
of  a  concentrated  charge  bigger  than  I  kg  when  the  scaled  charge  height  is 
about  4  ft/lb^/^  (1.6  m/kg^/3).  The  overpressure  produced  by  such  a 
charge  can  be  seen  by  referring  to  Figures  12,  46,  and  56;  assuming  a 
charge  weight  greater  than  1  kg  (2  lb),  overpressures  greater  than  10  Ibf/ 

•y 

in.  (70  kPa)  are  likely  to  be  applied  over  a  radius  significantly  greater 
than  the  thickness  of  the  snowpack.  However,  this  does  not  establish  a 
lower  limit  for  air  blast  release. 

Sonic  boom  studies  showed  that  avalanches  can  be  triggered  by  wide¬ 
spread  peak  pressure  of  2  Ibf/ft^  (0.014  Ibf/in. or  0.1  Pa).  This  is  the 
pressure  level  at  which  sonic  booms  begin  to  have  the  ability  to  crack 
plaster  and  glass  in  surface  structures.  It  is  equivalent  to  the  addition 
of  4  to  5  in.  (100-130  mm)  of  low  density  snow  (0.1  Mg/ra  )  on  slopes  of 
typical  angle  for  avalanche  release  zones.  In  air  blast  terras,  it  is  the 
pressure  that  might  be  felt  at  a  radius  of  40  m/kg^/^  (100  ft/lb^/^) 
when  a  charge  is  fired  just  above  the  surface  (say  0-5  ft/lb^/^,  or  0-2 
m/kg  1/3)  ^ 

For  practical  purposes,  it  has  been  suggested  (Mellor  1973)  that  the 

y 

target  zone  should  be  subjected  to  overpressures  of  at  least  0.5  Ibf /in. 
(3.5  Pa),  perhaps  by  firing  at  a  height  of  5  ft/lbl/3  (2  m/kg  1/3).  For 
a  charge  of  minimum  size  (1  lb,  or  0.5  kg)  this  gives  an  effective  target 
zone  of  diameter  100  ft,  or  30  ra. 
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